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shallow Phosphorous diffusion into p-type 
silicon was done, resulting in junction situated 
close to the surface. Diffusions were done at 
temperatures ranging from 850® G to lOGO^G. Junction 
depths varied from 0,5 micron to 2 microns. 

Impurity profile for each diffusion Wros evaluated 
employing anodic oxidation of silicon. The validity 
of movingfebundary solution of Fick’‘s second law was 
verified and the diffusion constants, of phosphorous 
were calculated. 
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IIST OP SYMBOLS 


a(x, t.) - Total Phosphorous Oonoentration at distance 

X from the surface for a diffusion time. t. 


0_^(x, t)&- 
GbCx, t) 


Concentration of the species A and B at x 
and for a diffusion time t. 


Os 

ObC^Co) 

and D-g— 

e 

K 

o< 

<r 

p 

PB 


Surface concentration of phosphorous 

Goncentration of the B species at the phase 
boundary Xq* 

Diffusion constants of A and ' B species of 
phosphorous . 

electronic charge. 

0B(xo)/Ga 

rate constant defined by xn . 

Bulk conductivity of the sample 
Bulk resistivity of the sample 
Sheet conductance of the sample 
Sheet resistance of the sa,mple 



INTRODUCTION 


Diffusion of acceptor and donor type impurities . 
into silicon is a "basic step in the fabrication of integrateff 
circuits. -A monolithic circuit can be arrived at by following 
any one of the three basic proces'^es (i) epitaxial -diffused 
process (ii) diffused-collector process and (iii) the triple 
diffusion process and in all these three, no matter what be 
the particular process followed, multiple diffusion of 
impurities is inevitable. As the performance of the indivi- . 
dual components in the chip and the isolation achieved by the 
diffusion when it is used for that (for eg in epitaxial diffu- 
sion pyocess and triple diffusion process) are of extreme 
importance in the final product^, the process of diffusion 
assumes -aii all important role in device fabrication. 

Once impurities are diffused into silicon wafer, 
to form iunctions, immediate thing one tries to find out is, 
how these impurities are distributed inside. This important 
Investigation serves two purposes. 

If one can characterise the distribution, as 
following a particular law during diffusion, any other time, 
given some known or assumed conditions the unknown parameters 
#f the system can be predicted. For example, if the distri- 
bution excatly Is error function like, for some required 
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junction depth one can find out how much time the source of 
impurities must be on, for other known parameters like tempe- 
rature background doping diffusion constant etc., So this 
necessitates the mathematical formulation of the problem, 
which gives out as a solution the expected distribution of 
impurities under the given conditions. Obviously this theo- 
ritical formulation must be experimently proved to arrive at 
correct predictions and so the necessity of evaluating the 
profile of impurities and comparing that with theoritical 
values, point by point, arises. This will be one aspect of 
the investigation. 

The other one is of design consideration. 
Knowledge of the concentrations and diffusion coefficients 
aid in the design of devices which will have prescribed 
electrical parameters. This aspect has been amply described 
in an article by Tannenbaum^, Device parameters like base 
width, base resistance, emitter efficiency etc. depend on 
diffusion conditions and these can be predicted closely 
well if the concentration, nature of distribution, and 
diffusion constants of the impurities are known beforehand. 

S'o in either way, knowledge of the diffusion, 
profile and the related diffusion constants of the impurities 
is essential and in this work an attempt is made to evaluate 
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the profile of phosphorous distribution in p type silicon 
under different temperature conditions, Some of the anamo- 
lies, "Which have been reported in the literatures in case of 
diffusion under some specific conditions, are of main concern 
in this -work. 
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CHAPTER - 2 


LAVJS OF DIFFUSION 


Th«re are two laws governing the process of 
diffusion in solids. 


First Law i J = 


Second Law 


6c 

5t 


- D 


b c 
b X 



dx^ 


The First Law governs the flux J and concentration 
gradient b 0/ dX where as the second one is concerned with 
the spatial distribution of impurities with respect to time. 
These two laws put together are called Fick*s laws. 


The law of interest to semiconductor technologist 
is the second one as the solution of this law for the given 
diffusion conditions ie time, temperature and concentration 
describes the impurity distribution in the solid. Solution 
•f this equations comes out to be different depending upon 
the boundary conditions and a look at ordinary diffusion 
conditions will be helpful in ariving at these boundary 
conditions . 


In most of the cases what we have is a diffusion 
furnace whose central zone temperature can be kept constant 
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within + l°Gj inside vjhich one can keep a quartz tube to 
house the wafers. The dopant mixed with the carrier gas will 
be let in constantly into the quartz tube^ the amount of 
gases being recorded by individual flow meters. If a liquid 
source like phosphorous-oxy chloride is being used oxygen 
also hps to be used along with the carrier so as to facili- 
tate the dissociation of phosphor ous-®xy chloride into 

P^O to form a local source near the wafer and also to 
d 5 

protect the silicon surface. Through out the diffusion the 
same rate of flow gases is maintained through the furnace. 
This arrangement ensures a constant supply of the dopant near 
the surface and this becomes an important bound, ary condition 
in the solution of Pick's second law. Because no impurities 
of the same conductivity type are present in the vrafer at 
the start of diffusion, we ^an see that at t=05 C= 0 . 

Under these conditions solution of second law 
can be seen to be^ 

G(x, t) = C3 erfc x/2v/Dt (l) 

where G(x, t) is the concentration of the diffused impurity 
at a distance x from the surface after a time t. Gg repre- 
sents the surface concentration; D and t represent the 
diffusion constant and time of diffusion respectively. 

Though this error function solution describes satisfactorily 
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the Infinite source diffusion in some cases, there are some 
cases where the distribution does not follow exactly this 
solution. It has been found that a moving boundary concept 
has to be introduced alongwith two diffusion constants into 
the boundary conditions to arrive at the diffusion profile 
under these special cases. These special conditions and 
how a moving boundary concept has taken shape are the subjects 
of the next chapter. However the solutions of the diffusion 
equation under thevce conditions are mentioned below. 

In this approach it is assumed that there are 
two species of phosphorous diffusing into silicon preceeded 
by a region where there is no change in the concentration 
ef impurities. A model fitting these assumptions is shown 
in Figure 1, 
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It consists of three distinct regions^ a constant 
surface concentration region followed by a transition region 
and, then the normal diffusion region. The assumed moving 
boundar ,7 is shown to be lying at a distance from the 
surface and for getting solutions of the dif'^'usion equations 
it will be assumed that the position of the phase boundary is 
directly proportional to the time the surface is exposed to 
diffusion} ie Xq =o<t. This boundary condition can be easily 
verified in practice and the validity of this is discussed 
by Tsai"^, The transition region consists of the two types of 
diffusing species A and B v/here as the normal diffusion 
region consists of only one species B. It is also assumed 
that for a given temperature the concentration of the fast 
diffusing species B is constant at the phase boundary and 
the ratio G3(xq)/Cs Is defined as Under these conditions 
one can write down the two diffusion equations for the two 
species. They are, 

O C^/c)t = ( 6 ^ 0 ^/ c'x^) 

and OCg/at = Dg (6^C-g/ dx^) 

where and Dg are the diffusion constants of 
the two species A and B. and Gg are the concentration 
of- the species at distance x from the surface. 
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Boundary conditions can be noted to be, 

Ga (x, O) =0 
Cg (z, o) = 0 
C_^ Coo jt) = 0 

Cg (® jt) = 0 

Cg (xq) = ri.Gg and C^^(Xq) = (l-lOGg 
Solutions to the tuo diffusion equations with the above 

4 

mentioned boundary conditions ares 


G^^Cz5t) = 


■. a .T-Kl 

2 


Cg exp 


(x- -'at )j\erf c S±.9^.t + erfc-?-i^.^~ 


2D-. 




2D, 


and 


Cp(x,t)= “ Cg exv\ -^^f-Cx-oct) 

® ® LyoB 


1 eric 

4- o T* *P r' 

- 8t-Xtl 

i L. 2 Dgt 

2 

Dgt J 


for a total solution 

for X < x^, ^ 

and X ^ x^ C = C^^Cx, t) + Cg (x, t) . 
This is the solution of Pick's law under the above 


mentioned conditions 


{ -P 
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CH.^TER - 3 

DEEP END 3 HALLOW DIFi^^USIOE? 

The lav.'s of diffusion and how far do they desc- 
rihe the practical diffusion profiles forms an important 
aspect of diffusion study, as lias already been Indicated. 

In the history of diffusion in semiconductors, 
first attempts were made by Fuller, Dunlap and others to 
get p-n junctions and they have reported the diffusion 
constants of some impurities of interest in Germanium 
around 1952^?® 

Dunlap diffused antimony in germanium and eval- 
uated the diffusion cons tents by tiro methods, one being radio 
active tracer method and the other p-n junction method using 
a thermal probe to detect the junction. 

In the same year Fuller from Bell Telephone 
Laboratories reported his resulfs of diffusion of ijs , Sb- 
In and Zn in Germanium follov/in''. a different method. Diffusion 
into two specimens of two different resistivities and ,^2 
was done for the same time under identical conditions. 

If the experiment is started v/ith p type speci- 
mens having gallium previously distrib’ited through out uni- 
formly, and i 7 ' 2 _ and cf2 conductivities of the two 
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specimens respectively^ then, 

rn = e p^ u and ^2 ~ e P2 

from this p, = and = 

" -L e .u 

Issumlny complete ionisation, 

C ,, = - — ™ and G „ 
e n ^"2 


62 

e 


0>2 
e M 


Diffusion into the'^e specimens was done and the 
depths and X2 upto which the impurity penetrates were found* 
Assuming that the diffusion has followed er^'or function so- 
lution, one must have 

C(x, ,t) = C, erfc and C(x2 > t)= C erfc — 

° 2 jDt 


But by the definition of .iunctioh these donor 
concentrations must be precisely same as the acceptor concen- 
trations existing there. 
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.-Cq erfc 
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= Cq erfc 


2 /bt 


Dividing; 




erfc 

2NiDt 


erfc 




2 'bt 

V 


0-, 

G'2 



11 


In this equation, the only unknown D can be eva- 
luated. These calculations are based on the assumption that 
Fick^s law was obeyed during diffusion. Attempts made to 
check this validity experimentally are given below. 

This law requires that when times required to 
reach a fixed arbitrrcry concentration at successive distances 
are measured, the distances must be proportional to the 
square root of the time. 

In the course of previous diffusions, it is 
evident that irrespective of the depth of location of p-n 
junction, the concentration of impurities near this layer 
will be alwa;'’'s the sarao. The only thing that is happening is 
that with prolonged heating this layer is going deeper and 
deeper into the substrate. This shift of the diffusion profile 
with time is well represented by Dunlap. 

It is possible to find this layer at a distance x^ 

^1 

if the diffusion is allowed for a time tp and at a different 
depth for a time t2» What all the above law requiret 



Figure 2 
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is that X.- ti and t.'x:. to or if we draw a graph between 
Ol" - 22 ~ 

O 

X and t, it must be a straight line . 

This is seen to be so from Figure 2. &‘o it is 

evident that diffusion .law is follovxed closely and that one 
can characterise the diffusion under these boundary conditions 
by the solution 'obtained from Fick’s law. 

Samilar attempts were made later by Fuller and 
others'^, to find oi.it the diffusion constants of different 
impurity elements in silicon. The assumption that the concen- 
tration of the impurity near the surface remains constant has 
also been checked up. 

In all these investigations and also in further 
investigations by various authors it was established that 
constant source diffusion can be described by erfc distribu*- 
tion. Experimental results of impurity concentration versus 
distance have been fitted into error function curves and the 
results were also satisfactory. 

Exposure of Anamolies: 

If one looks back into the graph for the veri- 
fication of Fick's law the diffusion times encountered were 
hundreds of minutes and the diffusion depths also were 
correspondingly high. But with the advent of Integrated 
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Circuit Technology. short diffusion times and so correspond- 
ingly small diffusion depths have become common. Ncw-a-days 
shallow phosphorous diffusion is normally employed for the 
emitter of an n-p-n transistor. 1 typical value of the depth 
of this emitter diffusion will be around one micron®. 

Follox^ring the vrork of Fuller and Ditzenberger , 
Tannenbaum continued the investigation of diffusion of 
phosphorous in silicon and reported that for relatively short 
diffusion times, the diffusion profile deviates considerably 
from the simple error function solution. It was found that 
electrically activa phosphorous concentration remains constant 
for a few thousand angstroms from the surfacelO. 

At the same time, a real check on the validity 
of error function distribution during constant source diffusion 
was made by Subashiev and others at the Academy of Sciences 

In this article they have suggested the existence 
of two diffusion constants instead qf one and they have 
conclusive!}’" proved that for short diffusion tiiaet the 
impurity profile can not be described by error function 
solution. 

In the present investigation, an attempt to 
check theoretical formulation for the experimentally observed 
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distribution of shallow phosphorous diffusion and the moving 
boundary concept is made at different temperatures. Practical 
diffusion profiles were obtained. Diffusion constants at 
different temperatures were evaluated from these curves and 
on the basis of the solutions given in cahpter 2, the concen- 
tration of impurities with respect to depth was recalculated. 
This calculated concentration is plotted back again on the 
same sheet} to compare it with the thecritical profile. 
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CH.iPTSR - 4 
EXP5RIMSNT/.L WORK 

D oping S y stem : 

The doping system vras set up in the top most 
shell of the teal-cwood almirah specially got made for this 
purpose. A plan of flovr of gases and the other accessories 
is given in Figures 3. 

Presently three flow meters are being used two 
for nitrogen and one for oxygen. Cne of the two flow meters 
for nitrogen, the larger capacity one, is used for measuring 
the flovr of nitrogen used as a carrier gas and this can be 
maintained at any required flow rata having a. maximum value 
cf 1450 co/min, for sapphire ball. The stainless steel ball 
can extend the range of the flow meter abox'-e this, but in , 
normal diffusion runs^l450 cc/min of carrier gas in quite 
sufficient. 

The other flow meter for nitrogen cf smaller 
capacity, .which is connected to the solonoid value is meant for 
bubbling through the liquid dopant, which is POCI 3 in this 
case. There are two possible states for the solenoid valve, 
normally closed (NC) and normally open (NO). When there 
is no power input to the valve NO side will be kept open 
for the gases. When power is on^ functions of these gates 
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are .just reversed. NC turns open and NO becomes shut. 

Flov of oxygen is necessary during diffusion 
both for helping the dissociation of POGlg into P 2055 so 
as to form a local source near silicon surface, from which 
actual phosphorous dif:^'usion takes place and also for pro- 
tecting the silicon surface. 

Liquid dopant is kept in the bubbler shown in 
the diagram. When diffusion is to take place, the solenoid 
Valve is kept on so that the side of the solenoid valve 
connected to the bubbler is open. Autom-i.tically, the other 
side of the valve which is connected to the main gas stream 
will be kept shut, so that no additional nitrogen can pass 
through the bubbler. This nitrogen which is bubbled through 
the dopant takes along v/ith it thi dopant and gets mixed with 
the main nitrogen and oxygen flows n.ar the entr.ance of the 
furnace tube. 

Temperature of the furnace is brought to the 
required value by setting the controller properly and it is 
advisible to measure the actual temperature inside the 
furnace with a platinum + platinum, 105^ Bhodium thermo couple. 
It is observed that there is not much variation in the set 
temperature with moderate changes of carrier gas flow. 
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Wien the furnace is brought to the required 
temperature the two nitrogon flow meters are adjusted to the 
required flow rates. Solenoid valve is kept off during this 
time. Cleaned wafers are then loaded carefully into the 
wafer-carrier and it is pushed approximately to the centre 
of the furnace tube. Now the oxygen supnly is turned on 
alongwith the solenoid valve, fliis is the time at which 
diffusion starts. 

After allowing the nitrogen to pass through the 
bubbler for the required time of diffusion, solenoid valve 
is put off, thus cutting out the supply of dopant. Oxygen 
flow can be continued for a couple of minutes more to grow a 
thin oxide layer over the diffused, surface which v/ill offer 
protection from the surroundings. Then both main nitrogen 
and oxygen flows can be cut off and the wafer unloaded. This 
can be stored either in methanol or in a clean^ covered 
petri-dish until further investigation is started. 

Clean ing Pr ocedu r e s Prior to Di ffusion 
For Glass-ware: 

All the beakers, petri-dishes etc are first 
washed with tap water throughly and then kept in hot EOH 
solution for about five minutes. After this they are taken 
out and washed in tap water and then- kept in chromic acid 
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solution for about five minutes. Then they are taken out, 
washed in tap water, DI water and finally st@']?ed in a decs- 
cator unit further use. 

Gleaning of silicon wafers : 

The 1 ohm-cm p type wafers used in this investi- 
gations were chemically polished and mirror finished, so 
that they can be directly effected to any one of the follow- 
ing cleaning procedures, 

Ist Method; 

Wafers are first cleaned ultra sonically in 
trichloroethylene for about five m5.nutes. 

Heating of those wafers is done successively in 
trichloroethylene and in a solution of equal parts, by 
volume, of H 2 O 2 , NH 4 OH and DI water for 30 minutes. 

Those are then taken out, rinsed in DI water 
and immersed in aquaregia for one minute* 

xifter washing them in DI water, wafers are 
immersed in HF for 30 seconds. This is followed by a rinse 
in. DI water. The two operations indicated in the above 
paragraph and this are repeated five times before a final 
through wash is given in DI water. Those wafers are then 
stored in methanol. 
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Second Method: 

Wafers are first cleaned tiltrasonically in 
trichloroethylene, heated for five minutes in it and then 
■washed throughly with DI water. 

1 solution of DI water, Hydrogen peroxide and 
HH 4 OH in tho ratio of 7:2:1, hy volnme, is prepared and heated 
to a temperature of 80®G, V/afers are then put into this hot 
solution for ton to fifteen minutes. After this, they are 
transfer ed into another beejeer and washed throughly in DI 
water. 

Another solution consisting of DI water, H 2 O 2 
and HCl (35^) is prepared in the same ratio of 7:2:1, by 
volume. Wafers are put into this solution, when it is heated 
to a temperature of 80°C. After keeping them for ten to 
fifteen minutes, they are taken out, washed throughly in 
DI water and stored in methanol. 

Of the two procedures tho second one was found 
to be less time consuming and most of the time, this was 
followed. 

Evaluation Procedure; 

Various methods have been employed to get the 

concentration of impurities as a function of diffusion depth 

from 

■ into the silicon wafer, /the surface. In any method, it is 
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necessary to remove thin layers and observe the electrical 
or other physical parameters after each removal of layer and 
then relate these observed characteristics to the concentration 
of impurities there. V'/hen the diffusions were deep^ this 
removal of layers was done by lapping with carefully designed 

jigsl2. 


Chemical etching can also be used for removing 
thin layers but here the problem is to control the thickness 
of layer removed each time. 

Anodic oxidation followed by dissolving the 

^?i02 layer in HF offers a very precise control of the 

thickness of the layers removed with very good repeatability, 

Specially for shallow diffusions this is much more attractive 

because one can remove very thin layers of the order of 
o 

100 A . By varying the voltage applied to the anodic cell, 
thickness of the layer removed, can be varied. As in 
this investigation shallovj diffusions are taken up this 
method of removal is employed s^nd is found t^ be most convenient 
and time saving. 

For this purpose, a special equipment was 
designed and is shown ta. figuir® 4, diagramatically. It 
consists of stainless steel piece with a provision to 
connect to the vaccuum pump on one side and on the other side 
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a graphite piece. can be fitted into this with the screw 
shown, w'^ilicen piece can be held "fight to this graphite 
piece. S'tainlo-ss steel bar is connected to the positive 
of a DG supply so that in the process of anodisation the 
silicon piece hecora. s positive and acts as the anode.- 

The catho-de of the call is made up of a 
platinum wire, wound in tho form of a circular spring. 

Even though, at present this is giving almost uniform 
oxide layers, employment of a platinum foil would be more 
appropriate. This would help in establishing a more 
uniform ionic current density all along the wafer surface. 

Both tho anode and cathode are fitted into a 
horizontal bar (with appropriate isolation), permanently 
fixed in the ring of a travelling microscope. By loosening 
the screw of the vertical slider the anode-cathode assembly 
can be moved up and down and fixed at whatever height It is 
needed. An inherent advantage of this fixed anode-cathode 
arrangement is that, the height of the. electrolytic column 
lying between them is same all along the exporiment'^tion and 
any variations in ionic current due to variation of electro- 
lyte resistance, from one setting to another, -are minimised. 

A 0.04 normal solution of ethylone glycol and 
Potassium nitrate has been employed as the electrolytic 
solution^^y with a small amount of water. 



22 


Electrolytic solution with the above composition 
is taken in a beaker and kept on the pl'^.tform of the travel- 
ling microscope. Anode is connected to the positive of 
a 500V, 100mA power supply and the cathode is connected 
to the negative of it. Vaccuum end of the stainless-steel 
piece is connected to a vaccuum pump, through a trap. The 
silicon piece to be oxidised is held under the graphite 
piece by the vaccuum pump. Electrical ccnnections are shewn in fi 

There is. a DPDT switch in the circuit one side 
of which is crinnectod to the anode-cathode assembly and 
other side is connected to resistivity measurement circuit. 
Keeping the switch in anodisation position, voltage is set 
around lOOV, For all the anodisations voltage between anode 
and cathode is maintained at 100+2 volts. After setting 
the power supply voltage to this value, the switch is 
thrown to the otherslde so that there is no supply to anode 
and cathode. Now the anode-cathode assembly is lowered such 
that only the bottom surface of the silicon piece comes into 
contact with the electrolyte. When this is done the slider 
is moved just a little bit up so that the liquid forms a 
meniscus with the silicon piece. This will lessen the 
possibility of any liquid creeping to the other side of 
silicon piece and thus gradually establish a direct 
connection between the electrolyte and the graphite piece. 
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This aspect of ad.lustment must he given sufficient importance 
because once direct conduction starts between electrolyte and 
graphite, the current through the circuit will not decrease 
with time as expected, and silicon dioxide layer ceases to 
form on the silicon piece. 

Once the silicon piece is adjusted as above, 
switch is thrown to anodisation side. Initial current 
through the electrolyte will be around 50mA, for the condi- 
tions in this arrangement, and this current gradually starts 
falling off due to the formcotion of silicon dioxide on the 
surface of the silicon piece and the corresponding increase 
in path -re si stance. In these experiments anodisation is 
stopped when the current falls off to a small value like- 
5 mA. 

Yaccuum pump is’ put off and the silicon piece 
is removed from the graphite piece with the help of 
tweezers. One must be careful here^ because, if a little 
bit of vaccuum is still left in the lines and the silicon 
wafer is tried to be taken away forcibly, it may break 
into pieces. That is why it will be better to have a 
stand-by sample diffused under the same conditions as the 
first one and should any thing happen to the first sample, 
evaluation can bo started on the second one, thus saving • 
a lot of time needed for cleaning and diffusing a fresh 
piece again. 
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. Tlie anodised piece, after removing from the 
graphite piece is ■well 'v/ashod in tap ’Water and then dropped 
in concentrated HP (48^) for about 30 seconds and again 
"Washed in water. This disealves the silicon dioxide layer 
formed on silicon due to anodisation and exposes a fresh 
layer of silicon on which 4-polnt probe resistivity 
measurements can be made. 

Determination of the thickness of silicon layer removed s 

A rough ostlm?.to of the thickness of the silicon 
dioxide layer formed can be got from the colour of tho 
film. After washing the anodised piece in water, it is 
allowed to dry and the colour can be seen under vertical 
flouracent light. For the conditions employ ■:rd in these 
experiments ie, lOOV, 50nA to 5nA, an^^ for the particular 
anodising solution used, colour of the film is always 
observed to be Brownish. This colour is also mostly uniform 
indicating that, the film of tf5i02 is uniform. By refering 
to the colour charts^^, one can see that this corresponds 
to a film thickness of 500 As it is known that for the 
formation of 100 thick Pi02 layer only 44 n of silicon 
Is made use of, this 500 2 . thick 5‘i02 layer corresponds 
■to 220 i of actual silicon removed^^. 
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This is only a rough estimato of the silicon 
layer renovecl and a more accurate method would he to find 
out hy thei weight difference. For this, the weight of the 
sample before any anodisationj is taken exact upto six decimal 
places using a microbalonce . After a knoi^m number of 
oxidations are done, the weight is taken again. The diffe- 
rence in weight gives the weight of total silicon removed. 

As each anodisation is done under exactly the same conditions., 
this weight difference by the number of anodisations gives the 
Weight of each silicon layer romoved.. As the area of the piece 
will be known, assuming the density of silicon to be 2.3 gms/cm^, 
thickness of e'^.ch silicon layer rcraoved can bo computed. 


In the Investigations a 287 mils x 287 mils 
square sample was used, so that, if jCox is the thickness of 
each layer- removed. 


( 287 ) 2 ( 25 . 4 ) 2 ( 10 - 4)2 


.•Nx X No of anodisations= 


Weight difference 
Density of silicon 


Shieet-Reslstance measurements : - 

After taking out the diffused sample from the 
furnace it is washed in HF for 30 seconds to remove the 
silicon dioxide formed on the top of it and first sheet 
resistance measurements are made on this surface. 
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For this 5 the DPDT switch is kept in Resistivity-* 
measurement position in the circuit. As surface sheet 
resistance will be very lessj of the order of a few ohrns/sq., 
current in the order of 20-50 rnAj through the current probes 
will be convenient to establish a potential difference 
accurately noticeble by the potentiometer, connected across 
the centre two V'-ltasc probes. The series resistance in the 
circuit and the volto.ge out put of the power supply are 
adjusted in such a way that the required current passes 
through the current probes. The potontionetar is standardised 
and the induced voltage betiNrean the two voltago probes is 
noted. Now the current is changed to some other value, by 
adjusting the power supply voltage and again the potential ■ 
difference between the voltage probes is noted. The Y/.^ 
ratio obtained with these two current levels must be approxi- 
mately same or otherwise, readings at different current 
levels must be taken such that for two consequtive current 
levels the 7/j ratio turns out to be approximately same. 

Another point to be carefully noted in this 
context is the position nf the sample under the probes. 
Variation in probe position on the sample from one set of 
readings to another set may make difference. To avoid this, 
the perspex cube on which the sample is being mounted is 
always kept exactly between two perpendicular guide lines 
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and then the silicon piece is adjusted on the top of this 
cnhe such that the same edge of the silicon piece is in contact 
with the fron edge of the perspex cuhej as shown in Fig.5 B 


Figures 5B 

Point A is marked on the pers-^pex cube so that the lefthand 
side corner of the front edge of the silicon piece, is always 
set at this point. Finally tho whole adjustment must he 
made in such away that the four probes of the header lie on 
straight line parallel to the edge and passing through centre 
of the top face of silicon piece, 

The process of finding out the sheet resistance 
is followed by anodisation ^washing and dissolving the Si02 in 
HF and again sheet resistance measurements on the fresh layer 
exposed. 'This process must be continued until the p-n 
junction is reached, which can be detected by a hot probe. 

Correction Factors 

Square samples having a side of 287 mils have 
been used for all the diffusions. Thickness of the piece 




is found to be 9 mils. As diffusion takes place from both 
the faces correction factor corresponding to this situation 
are taken into account. This factor comes to be 4.5414 for 
the samples under consideration. Multiplication of Y/j ratio 
at oach stage gives the sheet resistance at that point. 

Getting impurity concentration; 

When sheet conductance and distance into the 
sample are plotted on a linear graph, the curve is found to be 
linear upto a point. This corresponds to a constant 
concentration region vrhose bulk conductivity is given by 
d<rs/ gj X. 

For getting the bulk conductivity at other points 
procedure suggested by^^ Evans and. Donovan is followed. 

According to that, 

0.4343 fgCx) 

as ih^io ("s 

A curve has been fitted for each evaluation, ■ 
between logarithm of sheet resistance and distance. The 
derivative of that at all the points is compted which gives 
the denominator of the above expression. As Ps^x) at each 
point is known, the complete expression ’ ca.n be calculated, 
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which gives the bulk resistivity at each and every, point. 

Reference to Irvin curve is then made to 
find the impurity concentration of phosphorous corresponding 
to this resistivity. All this procedure is illustrated by 
an example. 

Getting diffusion constants i)^ and Dp:- 

Fron the curve got by plotting concentration 
versus distance two points arc taken one near the .Junction 
and the other one in the transition region. Near the junction^ 
as we have only one diffusing species, faster type, this 
concentration corresponds to that species. As the tine of 
diffusion and constant concentration region are known, the 
rate constant can be calculated. By extending the normal 
diffusion region curve backwards to the phase boundary 
CgCxo) can be found and thus K, which is G3(xo)/Gg. S?o in 
the equation, 

Gg(x,t)= - Gg exp ~ -ext) epfj + erfo . ~ — 

2 SDb SgDjt SpBtV 

every thing tScept Dg is known and so Dg can be calculated i 
For this a computer pro gr. am has been designed and is exp-- 
lained in the next chapter. 
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S-iiailarly, for gettinr. D, the concentration 
of the impurity at a point in the transition region is 
noted. This will be the total concentration of the two 
species present there. To get the concentration of A 
species only, th^ concentration of B species is found from 
the backward extended curve and deducted from total concen-^ 
tration. Knowing G^«^(x,t) as well o.s o.ll other parameters, 
is calculated. 

Gheckiiif; un theory with experiment 

To verify whether the above mentioned solutions 
exactly describe the nature of distribution, the computed . 
values of D/^ and Dg are back substituted into thoir original 
equations and concentrati..;ns G^pCx,t) and GgCxjt) are 
calculated at each point of interest. These two concen- 
trations are added to give the total concentration of impurity 
at each point and this is plotted again on the same graph* ; 

One can see that there is good agreement between the twe?, 
Besults are shewn graphically frem dPigures 6 t® 15, 
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MODEL WIJJJLTIO^ 
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A typical evaluation example at 900®G is given 
belo-w in detail. 

Wafer started with. : 1 ohn - cm p type 

Dimensions of the Silicon piece 287 x 287 x 9 mils. 

Diffusion temperature 900°C 

Diffusion time J 1 Hour. 

Main nitrogen gas flow; 1450 cc/mln 

Nitro.~,Gn flow through the bubbler: 30 cc/min 

Oxygen gas flow : 75 cc/min 

Initial weight ^^f the sample ) 

after diffusion and bef-.-re any ) 0.027865 gms 

anodisation is done ) 

Anodising voltage = 100 + 2 V 

Initial current through electrolyte = 50 + 4 mA 

Final current when anodisation is stopped = 5mA. 

Number of anodisations required) 
to reach p-n junction ) = 42 

Weight after 25 number 'of anodisations = 0,027803 gms 

o 

Thickness of silicon layer removed each time = 203 A 

As has already been mentioned, after each oxidation 
the sheet resistance was measured with Four point probe 
equipment. Bheet conductances corresponding to the first 
nine readings almost fall on a straight line, as shown in 
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the linear graph, and bulk conductivity of this surface layer 
is the slope of this straight line. This can be seen to be 
1/4. 2xlO”^a_ppPQ3^^ which corresponds to a concentration of 
3 X 10^®/ cm^ of phosphorous atoms. Bulk resistivity at all 
other points was found out using the expression already 
given. For calculating this, it is necessary to know the 
slope of logarlthra of sheet resistance at each point, and so 
a curve is fitted between x and pg with the help of computer 
and then the slope is computed. Those details are given 
below for 900** G diffusion. 


The following polynoninal was fitted to the 
Values of x and log (pg) 

LogioCpg) = ®*8Q8''i63^ + 4^3744 x 0w62744 x^-8. 9449x^-1.5, 957x^ 
where x is in microns 

This was differentiated at all the points to give the 
derivative and then RBO was calculated using the standard 
expression. Then from Irvin’s curves concentration 
corresponding to each value of RHO was noted. 

This concentration was plotted against x 
on a semi -log paper and is sh^vm as. the experimental 


curve 
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X 

(microns) 

(ohms/ sq 

^(logioPs) 

dx 

P 

(ohms -cm) 

G 

(atoms/cm®) 

0.20 

43.8 

37987.2 

5x10“^ 

2.5x1033 

0.22 

52.0 

36892.9 

6.1x10"^ 

1.8x102P’ 

0.24 

62.5 

35773.2 

7.6x10”^ 

1.2x1033 

0.26 

75.3 

34645.2 

9.4xl0~^ 

9x10^^ 

0.28 

87.5 

33526.0 

1.1x10“^ 

7x10^® 

0.30 

99.8 

32433.0 

1.3x10"^ 

6x10^3 

0.32 

116.5 

31383.2 

1.6x10“^ 

4.5x10^^ 

0.34 

135.0 

30394.0 

1 .9x10 

4x1 0^® 

0.36 

151.0 

29482.5 

2.2x10-3 

3.5x10^^ 

0.38 

171.0 

28665.8 

2.5x10“^ 

2.7x10^® 

0.40 

194.5 

27961.3 

3x10 

2.4x10^® 

0.42 

220.0 

27386.0 

3.5x10-3 

2x10^^ 

0.44 

250.0 

26957.3 

4.0x10-3 

l.Sxlol® 

0.46 

290.0 

26692.2 

4.7x10-3 

1.4x10^3 

0.48 

329.0 

26608.0 

5.4x10"^ 

l.lxlOl® 

0.50 

37S*C 

26721.9 

6.0x10-3 

1,0x10^® 

0.52 

434. « 

27051.1 

6.9x10“^ 

8.5x10^3 

0.54 

505.0 

27612.8 

7.9x10-3 

7x10^^ 

0.56 

538.0 

28424.2 

8,2x10-3 

6.5x10^3 

0.58 

605 . 0 

29502.5 

8.9x10-3 

6x10^3 

0..60 

715.0 

30864,8 

1,0x10“^ 

5.0x10^2 

0.62 

852.0 

32528.4 

l.lxl0~2 

4.0x10^^ 

0.64 

9f?0.0' 

34510 • 5 

1.2x10-2 

3.5x10^^ 


Gontd. on the next page. 
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• X 

(ni crons) 

fs 

ohns/sq 

^ (logio pc) 
dx 

1 

(ohms -cm) 

C 

(atoms/cm^' 

0.66 

1100.0 

36828 .3 

1.3x10"2 

3x10^^ 

0.68 

1300.0 

39498.9 

1.4x10-“2 

2.8x1018 

0.70 

1610.0 

42539.7 

1.6x10“2 

2.0x10^^ 

0.72 

1950.0 

45967.7 

1.8x10’*^ 

1.7x10^® 

0.74 

2600.0 

49800.1 

2 . 2x10-2 

1.3x10^® 

0.76 

3010.0 

54054.3 

2. '4x10 "2 

1.0x10^.® 

0.78 

4380.0 

58747.3 

3.2x10“^ 

6.0x10^'^ 

0.80 

5430.0 

63896.4 

3. 7x10 

4.5x10^'^ 

0.82 

7230.0 

69518.7 

4.5x10"2 

3x1 0^'^ 


From th 

e experimental 

curve two points were 

taken, one at x^ = 0.24 and another at X 2 = 0.70 microns. 

' Value of K - ^ = 0.24 

3.0 X 1020 

The constant surface c'^ncentration region 

extends upto Xq = 0.18 x lO""^ cm. 

Sk) 

, from Xq 

= cK t , ■ <x, = 

0.18 X 10"'^ 
3600 

=5 X 10’® 

/it 

Xl, C(xi) 

s 1.2 X 10^0 

and from the. ex; 

tended 


curve GbC^I^ = 5 x 10^® , so thot (x]_) - 7 x 10^®. 


This value o,f Ga (xi) and' the above given values of k, (x, j 


t 
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and X were sulDstltuted in- the equation for ^ to get 

the value of the unknown tern D.*» 

18 

Ag-ain, at X 2 j G3CX2) = 2 x 10 ^ and this value of 

concentration at X 2 and all the above values of k, oCand t 
were substituted into the equation for Cg(xj t) to get the 
value of Dg. 

The values of % and Dg calculated in the above 

manner are, 

D. = 1.27 X 10“^'^ 
and Dg = 6.2 x 10”^"^ 

These values were substituted back into the two 
original equations and the concentrations Gj^Cx, t) and 
Gg (x, t) were calculated at all the points. The sum of 
these two gives the total concentration CCx, t) . These 
are tabulated for this casey as shown cn thes.- next page. 

This concentration, calculated frorx the theori- 
tical equations, was plotted again with respect to depth and 
is shown as the theoritical curve, ks the experimental and 
theoritical curves are in good agreement, it can be conaluded 
that the' concept of two diffusion constants is valid., A 
discussion of the results is given in chapter 6 in a more 
detailed way. 


X 

li crons) 

atons/ cm'^ 


t, 

atoms /cm^ 

) 

C=C_^Cb 
atom’s/ cn^ 


0.20 

0.153832E 

21 

0.652550E 

20 

0,219087E 

21 

0.22 

0.103781E 

21 

0.591042E 

20 

0.162885E 

21 

. 0.24 

0.700000E 

20 

0.534811E 

20 

0,123481E 

21 

0.26 

0. 471956 E 

20 

0.-183291E 

20 

0.955247E 

20 

0.28 

0. 317975 E 

20 

0.436005E 

20 

0.753981E 

20 

0,30 

0.213974E 

20 

0.392552E 

20 

0.606527E 

20 

0.32 

0. 143711 E 

20 

0.352595E 

20 

0.496306E 

20 

0.34 ■ 

0.962369E 

19 

0.315848E 

20 

0.412085E 

20 

0.36 

0.641692E 

19 

0 . 282073E 

20 

0.346242E 

20 

0.38 

0.425309E 

19 

0e251066E 

20 

0.293597E 

20 

■ 0.40 

0.279631E 

19 

0.222e50E 

20 

0.250613E 

20 

0.42 

0.181952E 

19 

0.196671E 

20 

0.214866E 

20 

0,44 

0.116872S 

19 

0.172988E 

20 

0. 18467 5E 

20 

0.46 

0.739074E 

18 

0.151474E 

20 

0.158864E 

20 

0.48 

0 o458906E 

■18 

0.132005E 

20 

0.136594E 

20 

0.50 

0;279048E 

18 

0.114465S 

20 

0.117256E 

20 

0,52 

0.165757E 

18 

0.987388E 

19 

0.100396E 

20 

0.54 

0.959627E 

17 

0.847110E 

19 

0.856707E 

19 

0.56 

0.540328E 

17 

0.722674E 

19 

0.728077E 

19 

0.58 

0.295338E 

17 

0.612933E 

19 

0.615886E 

19 

0.60 

0.156446S 

17 

0.516742E 

19 

0.518306E 

19 

0.62 

0.801962E 

16 

0.432964E 

19 

0.433766E 

19 

0.64 

0.397313E 

16 

0.360478E 

19 

0.360875E 

19 

0.66 

0.190027E 

16 

0.298188B 

19 

0.298378E 

19 

0.68' 

0.876567E 

15 

0.245034S 

19 

0.245121E 

19 

0.70 

0.389650E 

15 

0 . 200000B 

19 

0.200039E 

19 

0.72 

0.166790E 

15 

0.162125E 

19 

0. 162141 E 19 

0.74 

0.687059E 

14 

0.130507E 

19 

0*130514E 

19 

0.76 

0.272215E 

11 

0.104313E 

19 

0.104315E 

19 

0.78 

0.103684E 

1-1 

0. 82778 2E 

18 

0.827792E 18 

0.80 

0.379506E 

13 

0.652123E 

18 

0.652127E 

18 

0.82 

0.133435E 

13 

0.509965E 

18 

0.509966E 

18 
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CHAPTER - 6 

DI?Cn?3'J0N OF THE RESULTS' 

From the plots at different temperatures, it is 
evident that, the distribution of impurities can be characte- 
rised by three distinct regions. As has already been indicated^ 
these are a constant concentration region, a transition 
region, followed by the nornal diffusion region* 

. Presence of constant surface concentration region 

has been established by many authors. Tannenbaum was the 
first to give a detailed account of this for phosphorous 
diffusion is silicon. This was followed by many reports about 
the existence of a constant surface concentration during 
Boron diffusion as well as under different conditions of 
phosphorous diffusion. 

There is always a difference between the concen- 
trations of phosphorous atoms, as measured by electrical 
means, and the total phosphorous concentration as can be 
obtained by radio tracer measurements. How a flat region 
arises near the surface can be explained in the following 
way . 

At high concentrations of phosphorous, diffusion 
introduces dislocation networks, poaching such densities 



as 5 X lO^cn/cn^ nerr th:.- surface. Then networks can be 
found to be roachino' such depths as 1.5 ji fron the surface^^. 
At still Croat ,jr distanc 'S these netw-^rks degenerate into 
individual dlsloca.tions and finally reaching the densities- 
present in the original sanple, before diffusion. This 
generation of dislocation network was found to start with 
a concentration of 10^^/cm^ of electrically active phospho- 
rous atoms. , These diffusion induced dislocations have 
predominant edge character and as -already there is a super- 
saturation of phosphorous atoms they precipitate along the 
dislocations and at other strain centres. That this is so, 
has b-oen experimentally verified by using electron micros- 
copy-*-'^. 

The fact that Impurity concentration in the 
range of 10^^/cn^ gives rise to dislocation formation can 
be seen in the follovring way. 

Because of diffusi.:n of impurities there is a 
mismatch between solute and solvent and the strain due to 
this mismatch will be remove'" by the formation of dislocations 
If tr denotes the stress induced due to this impurity 
concentration, for shallow diffusions, it has been shown 
by Thai^^, that this is related to the impurity concen- 
tration by the relation 
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where /3 
E 
h) 

and C = 


= P’olute lattice contraction coefficient 
= Young Modulus of silicon 
= Poisson's ratio 

concentrati'-ai of impurity 


If this induced stress exceeds the flow stress 

i 

of the naterialj given by hN® yu/w, 

whore, h is the magnitude o-f Burgr's vector 
N is the dislocation density 
j:u is the shear modulus 

and w is a numerical factor ranging from 3 to 5 
IT, in itself is related to the concentration of 
impurities by the simple relation !T = OpC where Q is a constant* 
So flow stress = b (Q fo ja / w. 

Plastic flow occurs when induced stress' exceeds 
this, ie. when. 

c >. JL . r »b (1 - | 2 

fi . wE -J 

Substituting the vo.lues of various terms in 
this expression, it can be seen that this concentration 
comes , to be about 5 x 10^'^/cn^. 


This value of concentration, where plastic 
flow begins ie dislocation networks begin to form, is in 
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■whore yS = F'olute lattice contraction coefficient 
E = Young Modulus of silicon 
'9 = Poisson's ratio 
and C = concentrati'-^n of impurity 

If this induced stress exceeds the flovj stress 
of the naterialj given by bN® yu/Wj 

whore, b is the magnitude o.f Burgr's vector 
N is the dislocation density 
p. is tho shear modulus 

o.nd w ir a. numerical factor ra.nging from 3 to 5 
N, in itself is related to the concentration of 
impurities by the sinplo relation !T = Q|3C where Q is a constant* 
So flow stress = b(QjQG)^ U / "w* 

Plastic flow occurs when induced stress* exceeds 
this, ie. when, 

G > . r .Uh (1 - >~ ) ' ~j 2 

S'ubstituting tho values of various terms in • 

this expression, it can be seen that this concentration 

no o 

comes , to be about 5 x 10-^ ■'/cn'^. 

This value of concentration, where plastic 
flow begins ie dislocation networks begin to form, is in 


confirnation with the c-oncentratioa value where the electri- 
cal rioasurements deviate from radio tracer measurements. 

Once the concentration of impurities exceeds this, the 
silicon matrix into which these impurities are introduced 
starts adiustinp itself by thermally-activated self diffusion. 
That is why the solute atoms will propagate to the dislocations 
and to ether discrete strain centres, and this results in 
prccir-itation. ’ fo during diffusion, the strained lattice 
attenr-ts to maintain in solution an impurity content of not 
mo-o) tho.n the solubility limit corresponding to the thermo- 
dynarmlc stability. This only can be measured by electrical 
methods and is equivalent to the total ph'^sphorous in true 
solid solution. ' 

Knowing the density of dislocations with respect 
to distance into the wafer, one can calculate the amount of 
precipitation at each and every point and by subtracting this 
precipit.ation profile from the ideal error function profile, 
the true curve for the electrically active phosphorous will 
be obtained. This phenomenon leads to a constant surface 
concentration. 

In the present investigations also a constant 
surface concentration has been observed. When the measured 
sheet conductance is plotted against the distance, a curve 
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results in which the slope is cons tint up to a particular 
distance. That this is so, is evident from all the graphs 
shown, at different tenperatures . :.s the slope of sheet 

conductance versus distance gives the bulk conductivity, 
after finding the slope, reference is made to Irvines curve 
to get the surface concentration. Whan ever there is a 
snail deviation of sd mo points from the straight line, it 
is believed that it is due to the variation of probes 
position on the wafer from one reading to another. The 
general trend of, the distribution is clearly linear near 
the surface j thus establishing the fact that surface concen- 
trati'ui ranains constant, upto a distance. 

It was assumed, while g.-otting the solutions 
for diffusion, equati'^ns, that a phase boundary exists and 
the end of this constant surface conce.itratlon region 
corresponds to this. 

21 

Investigation by S'-chmldt et al revealed 
that the precipatlons are in the form of a compound of 
silicon and phosphorous and were verified to be silicon 
phosphide. Once the dislocation networks cease t^' exist 
at tho phase boundary some of the precipitated material 
may take up substitutional places, thus contributing to 
another diffusion coefficient. If it is assumed that 
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molecular silicon phosphide migrates without dissociation, 
its diffusion constant will be lower than elemental phospho- ' 
rous. This has also been verified in these investigations. 

For all the cases value of the constant is smaller than 
Dg. Thus it i'^ probable that in the region immediately 
following the phase boundary both the species are present 
and the total elontrically active phosphorous concentration 
will be the sum of the two individual concentrations. 4s 
the normal diffusion region is reached, the concentration 
of 1 species falls of' rapidly and this infact can be seen in the 
model evaluation givon. 
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APPENDIX 

Two computer program packages, ^BIJLKRf and 'DIFGOF^ 
have been used to compute the bulk resistance and the diffu- 
sion coefficient (and related calculations) respectively. 

A brief outline of the methods adopted in them- is given below. 

^BULKR^ accepts as input (with Runtime format 
facility) the number of observations of X made, the values of 
X and the corresponding values of pg. The values of 'log-LoPg 
are computed using the built in FORTRIN function ALOGIO(Z). 

The values of X and log^o Ps 9-re fed to a routine 
^POLY*. This routine obtains a good polynominal fit for the 
sets of values given. Given a number n, the values of x and 
y, a routine CURVE ^ obtains an n"^^ degree polynomial for 
the Yj_n terms of x. The method employed is the Regression 
coefficient method. The Routine 'CURVE* computes the n + 1 
coefficients of the degree polynomial; also it computes 
the values of Y at the given values of x using the polynomial 
and compares them with the given values of x and obtains the 
standard deviation, *POLY* repeatedly calls ‘CURirE* for 
various values of n and selects that polynomial which has 
minimum standard deviation. The limits of n between which 
all degrees are tried are controllable through appropriate 
input parameters. 
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^BULK* proceeds with the best polynomial obtained and 
computes the polynomial for the derivative of logio P s 
algebraic differentiation^ Using this the values of the 
derivative of log^ops are computed at various x. Then using 
the formula = 0.4343 the bulk resistivity 

is computed, A typical input to and output of the package 
'BULKR^ are given below. 

The input is as follows Format is 3 I 5 for the first 
card to give the number of observations the limits of the 
degrees of polynomials to be triedj lOE 8.2 for values of 
X and Pg. 


19 3 6 . 

.lOE-4, .12E-4, .14E-4, .16E-4, .18E-4, .20E-4, .22E-4, .24E-4, .26B-4, 
.28E-4, .30E-4, .32E-4, .34E-4, .36E-4, .40E-4, .42E-4, .44E-4 j ,46E-4, 

•136E+3, .172E+3, .208E+3, .246E+3, .280E+3, .330E+3, .377E+3, .554E-H:, 
.627E+3, .815E+3, .957E+3 , .131E+4, .272E+4, .754E+4, .122E+5, .192E+5. 


The output is as follows: 


THE FOLLOWING ARE Y/iLUES OF LOGiq 


,213354E + 1 
.251851E + 1 
.279727E + 1 
.328217E + 1 
.408458E + 1 


.223553E +1, .2lll806E-»-l,.239094E+l, .244716 E4^-l, 
.267634E 1 .263649E + 1 .274351E + 1 

.291116E + 1 ,298091E + 1 .311661E +1 

.343457S + 1 .364345E + 1 .387737E +1 

.438217E + 1 
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POLYNOMIAL FIT BEGINS 
DEGREE =3 S = 0.023032 

DEGREE = 4 S = 0.01846« 

DEGREE 5 S * 0.014035 

DEGREE =6 S = 0.015431 

THEREFORE SELECTED DEGREE IS = 5 FOR WHICH S 0.014035 
THE FOLLOWING ARE COEFFICIENTS OF POLYNOMIAJQ OF LOG RHOS 


Xo 




0.249768E + 1, -0.162668E + 2 0.205037E + 3 -0.9608E + 3 

X4 X5 

,0.^958B + 4 -0.154960E + 4 

THE FOLLOWING ARE COEFFICIENTS OF POLYNOMI/iL OF DERIVATr/S 
OF LOG EHOP 


-0.162668E + 6 

% 

0.815834E + 8 


0.410074E + 1 


-0.774802E + 8 


-.0,288224E + 8 


THE FOLLOWING ARB COMPUTED, VALUES OF RHO 
X RHOS DERIVE 


RHO 


This table was given for 900°C case, in the model 
evaluation; Chapter 5. 
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^DIFGOF^ is package which organises the computations 
of diffusion coefficient ♦..It makes use of aroutine ^GOMPUT*’. 
The values of Xi, XgjT, Ci^XT, GEXT, K, ALPH4, ,G?? are all 
available for 'GOMPUT'. In addition 'GOMPUT' has 4 parameters 
I,J, DA, and BF. DA is the value of 0/^ in the equation, 


2 

(1-K) 



G.a (x, t) 


2 - erf 




erf 


x-SoCt 

J 


oCCx 
2 Da 


¥hen 1=0, the routine makes use of the above 
equation and using GAXT if I = 1, or CBXT if I = 2, computes 
the L.H.S. and R,H*S. of the above equatiorj* ■ DF is set to 
the difference between the LBS' and RHS, 'When J = 1 the 
routine computes the values of dAXT (I = 1) or GBXT (I =2) 
and DF is set to that value, 

tDIFCOF* takes in the values of X,T,CAXT, GBXT, K, 
ALPHA, GS, RMi;S, MIB as input, MAX and MIN are the range 
of values between which DA lies. By repeatedly making use 
of the routine ^'COMPIJT* with appropriate control parameters, 
the DA is set to MIN, RMAX and (M;*X + MIW)/2 and the 
differences between the LHB and EBB of the above equation 
are computed. As' the value of li^-RHSi is monotonically 
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decreasing logarithmic search is made to find the Talue of 
DA which satisfies the equation, Ev^erytime the difference 
is computed at the middle point of EMAX and BMIN and if the 
difference is positive RM^IX is set to the average; if the 
difference is negative RMIN is set to the average. Then the 
procedure is repeated with the new values of RM/JC and RMIN, 

This continues until the difference between the LHS? and RIF 
is O.IE + 8. The value of DA at this point is given as the 
accurate value. Bimilary DB is computed. 

Using the values of DA and DE and the given values of X, 
the values of GAXT, GBBT, and GX are computed by using the 
routine ’GOMPUT* in the J - 1 mode. The results are tabulated. 

Typical input to and output of the package are shown 
below: - 

Input given Xp, X 2 , CiJCT, CBXT, T, K, .U^PHA, CS, RM;DC, RMIN 
and the values of X are 

.12E-4, ,40E-4, .25E+20, .30E+19, 3600.0, 0.1, 0.222E-8, 
0.30E^-21, 0.10E~9, O.lOE-16 
The output is as follows: 

GALGULATION OF DA: 

RMAX = O.lOOOOOE-9, RMIN = O.lOOOOOE-16, R = 0.500000E-10, 

DF = -0.9O6822E0O 

RMAX = 0.500000E-10, RMIN= O.lOOOOOE-16, R =0.250 OOOE~10, 

DF=-0 .906527E00 

RMAX = 0.187007E-14, miN= 0.187006E-14, R = 0.187007E-14, 

DT^-0.242144E-7 
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THE VALUE OF DA SATISFYING THE GIVEN EQUATION IS FOUND TO BE 
0.187007E - 14, 

Sdmilarly DB also can be found out. 

THE FOLLOWING ilRE COMPUTED V/ILUES' 

X C.DCT CBXT GX 


Tabulation of those values for 900“G diffusion has 
shown in the model evaluation, Chapter 5, 

jILL routines ARE WRITTEN IN FORTRfiN IV. 'BULKRt ig 
a 200 statement program and its average execution time is about 
2 seconds for a problem ^DIFCOF* is an 80 statement program 
and its average time of execution is 2 seconds. 
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